LOCAL RIGIDITY OF INVERSIVE DISTANCE CIRCLE 

PACKING 



REN GUO 



Abstract. A Euclidean (or hyperbolic) circle packing on a closed triangulated 
surface with prescribed inversive distance is locally determined by its cone 
angles. We prove this by applying a variational principle. 



1. Introduction 

1.1. Andreev-Thurston Theorem. In his work on constructing hyperbolic met- 
rics on 3-manifolds, W. Thurston ([31], Chapter 13) studied a Euclidean (or hyper- 
bolic) circle packing on a closed triangulated surface with prescribed intersection 
angles. Thurston's work generalizes Andreev's result of circle packing on a sphere 
[H [2] • The special case of tangent circle packing (every intersection angle is zero) 
on sphere was obtained by Koebe [T3] . 

Suppose (E, T) is a closed triangulated surface so that V, E, F are sets of all 
vertices, edges and triangles in T. We identify vertexes of T with indexes, edges 
of T with pairs of indexes and triangles of T with triples of indexes. This means 
V = {1,2,...\V\},E = {ij | i,j G V} and F = {Aijk | i,j,k G V"}. Fix a vector 
G R' B ' indexed by the set of edges E, such that Q i3 G [0, ^] for the each ij G E. 
This vector is call a weight on (E,T). 

A circle packing on a closed triangulated weighted surface (E,T, O) is a con- 
figuration {ci,i G V} of circles such that the intersection angle of Ci and Cj is 0^ 
which is the prescribed weight on the edge ij. To obtain a Euclidean (or hyperbolic) 
circle packing for (E, T, 0), we start with a radius vector r = (r\, r%, riyi) G 
which assigns each vertex i G V a positive number r^. A radius vector r produces 
a Euclidean (or hyperbolic) cone metric on the surface as follows. 

In Figure[T] consider a topological triangle Aijk. One can construct a Euclidean 
(or hyperbolic) triangle Aijp such that the edges ip,jp have lengths Vi,rj respec- 
tively and the angle at p is ir — ©„. Let kj be the length of the edge ij in the 
Euclidean (or hyperbolic) triangle Aijp which is a function of Ti,rj. In fact, in 
Euclidean geometry 

hj = J rf ~+rf+ iTiTj cos Oy 

or in hyperbolic geometry 

lij = cosh - 1 (cosh 7*j cosh r.j + cos y sinh sinh r j ) . 

Similarly, we obtain Ijk , lu ■ 

As observed by Thurston 5T], under the assumption that Oy G [0, 5] for each 
ij G E, lij, Ijk and Iki can be realized as edge lengths of a Euclidean (or hyperbolic) 
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Figure 1 . Andreev-Thurston circle packing 



triangle Aijk in F. Gluing all of these Euclidean (or hyperbolic) triangles in F 
produces a Euclidean (or hyperbolic) cone metric on (£, T) with possible cones 
at vertexes of T. On this surface with a cone metric, drawing a circle centered 
at vertex i with radius r\ for each vertex i G V, we obtain a circle packing with 
prescribed intersection angles. 

At each vertex i E V, the cone angle is the sum of all inner angles hav- 
ing vertex i. Thus cone angles (oi, 0,%, o,\v\) are functions of the radius vector 
r = (ti, T2, !"|yi). We try to find a radius vector when cone angles are given. 
Andreev-Thurston Theorem answer the question about the existence and unique- 
ness of solutions of radius vector when cone angles are given. 

Theorem 1 (Andreev-Thurston). For any closed triangulated weighted surface 
(£, T, 0) with Qij E [0, j] for each ij E E, a Euclidean (or hyperbolic) circle 
packing for (£,T, 0) is determined by its cone angles up to Euclidean similarity (or 
hyperbolic isometry). Furthermore, the set of all possible cone angles form an open 
convex polytope in Rl^L 

For a proof, see Thurston [3T] , Marden-Rodin [T7j , Colin de Verdiere [5j , He [13] , 
Chow-Luo [9], Stephenson [20] . 

1.2. Inversive distance. H. S. M. Coxster [7] introduced the notion of inversive 
distance to describe the relationship between two circles in a Mobius plane. The 
notion of inversive distance generalizes the notion of intersection angle of two circles. 

The inversive distance between two circles is independent of geometry, i.e., the 
inversive distance of two circles is the same under spherical, Euclidean and hyper- 
bolic metrics. This generalizes the fact that the intersection angle of two circles is 
the same under spherical, Euclidean and hyperbolic metrics. 

In the Euclidean plane, consider two circles c\,C2 with radii r\,r<i respectively. 
In the relevant case of this paper, we assume that does not contain Cj for {i, j} — 
{1, 2}. If the distance between their center is I, the inversive distance between c%, c-i 
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is given by the formula 

I(ci,c 2 ) = — i K 

2rir 2 

In hyperbolic plane, consider two circles Ci,c% with radii ri,r 2 respectively. We 
assume that c; does not contain Cj for {i,j} = {1, 2}. If the distance between their 
center is I, the inversive distance between ci, c 2 is given by the formula 

cosh I ~ cosh n cosh r2 



/(ci,c 2 ) 



sinh rj sinh r 2 




Figure 2. the stereographic projection 

These two formulas are related by a stereographic projection. In Figure [21 con- 
sider the upper half space model of the hyperbolic space H 3 . The hyperbolic plane 
H 2 is realized as a totally geodesic plane in H 3 , where it is a unit hemisphere perpen- 
dicular to the infinite boundary C. The stereographic projection is the projection 
staring from point (0,0, —1), which sends objects in H 2 into C. Let Ci,c 2 be two 
circles in H 2 with radii r%, r 2 and the distance between centers I. The stereographic 



projection produces two circles ci, c 2 on the Euclidean plane 
the distance between centers I. Then we have 



I(ci,c 2 ) 



cosh I — cosh n cosh r 2 V 



with radii r\ , r 2 and 



I(ci,c 2 ). 



sinh ri sinh r 2 2rir 2 

We say that the stereographic projection preserves the inversive distance. 

When I(ci,c 2 ) G ( — 1,0), the circles ci,c 2 intersect with an intersection angle 
arccos/(ci, c 2 ) S (fj^)- When I(ci,c 2 ) G [0,1), the circles ci,c 2 intersect with an 
intersection angle arccos/(ci, c 2 ) G (0, -|]. When i(ci,c 2 ) = 1, the circles ci,c 2 are 
tangent to each other. When J(ci, c 2 ) G (1, oo), the circles ci, c 2 are separated. 

1.3. Inversive distance circle packing. Motivated by the application to discrete 
conformal mappings, Bowers-Stephenson [5] introduced inversive distance circle 
packing which generalizes Andreev-Thurston's intersection angle circle packing. See 
Stephenson [5D] (page 331) and Bowers-Hurdal [3] for more information. 

We will use the same notation as in Andreev-Thurston's intersection angle circle 
packing. For a closed triangulated surface (E,T), fix a weight I G K' B ' with 
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Iij G [0, oo ) for each ij G E. In this paper, we will not consider the case of obtuse 
intersection angles. (Rivin [18] and Leibon [15] obtained results for certain circle 
patterns on surfaces allowing obtuse intersection angles.) An inversive distance 
circle packing for (S, T, I) is a configuration of circles {ci,i G V} such that the 
inversive distance of two circles Ci,Cj is the prescribed number /y for the edge ij. 
To obtain a Euclidean (or hyperbolic) circle packing for (S, T, I), we start with a 
radius vector r = (r±, r%, r\v\) € rLJq which assigns each vertex i G V a positive 
number r^. A radius vector r produces a Euclidean (or hyperbolic) cone metric on 
the surface as follows. 

Consider a topological triangle Aijk G F. Assign a number £y to the edge 
ij G E, where /y is the distance of two circles with radii j-j, rj and inversive distance 
Iij. Therefore /y is a function of rj, rj. In fact, in Euclidean geometry 



Similarly, we obtain Ijk , hi ■ 

To guarantee that lij, Ijk, hi satisfy the triangle inequality for each triangle 
Aijk G F, a vector r = (ri, r2, must satisfy a certain conditions. For 
details, see section 2.1 and 3.1. 

For a radius vector r, if lij, Ijk, hi can be realized as edge lengths of a Euclidean 
(or hyperbolic) triangle for each Aijk G F, gluing all of these Euclidean (or hyper- 
bolic) triangles in F produces a Euclidean (or hyperbolic) cone metric on (£,T) 
with possible cones at vertexes of T. On this surface with a cone metric, drawing 
a circle centered at vertex i with radius rj for each vertex i G V, we obtain a circle 
packing with prescribed inversive distance. 

The cone angles [at, a,<i, a\v\) are functions of the radius vector r = (n, ri, r\y\). 
We are interested in the solutions of radius vector when cone angles are given. 

Our main theorem is the local rigidity of inversive distance circle packing. 

Theorem 2. For any closed triangulated weighted surface (E, T, I) with Iij G [0, oo) 
for each ij G E, a Euclidean (or hyperbolic) inversive distance circle packing for 
(S,T, 7) is locally determined by its cone angle up to Euclidean similarity (or hy- 
perbolic isometry). 

The local rigidity of Andreev- Thurston Theorem is generalized to the case of 
inversive distance circle packing. In our case, we are not able to obtain either the 
global rigidity or the description of the space of cone angles as in Andreev-Thurston 
Theorem. For discussion, see section 2.1 and 3.1. 

In section 2, we prove Theorem [2] in Euclidean geometry. In section 3, we prove 
Theorem [5] in hyperbolic geometry. 

1.4. Variational principle. Theorem[5]is proved by applying a variational princi- 
ple. The energy function of this variational principle is derived from the derivative 
of the cosine law of a Euclidean (or hyperbolic) triangle. Andreev and Thurston ob- 
tained their results without using variational principles. The variational approach 
to circle packing was first introduced by Colin de Verdiere [5]. Since then, many 
works on variational principles on circle packing or circle pattern have appeared. 




or in hyperbolic geometry 



lij = cosh 1 (cosh coshrj + -fy sinhrj sinhrj). 
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For example, see Bragger [3], Rivin [18], Leibon [15], Chow-Luo [9], Bobenko- 
Springborn [6J, Guo [IT] . Springborn [19], Luo [16], Guo-Luo [12] and others. 

Applying a variational principle, Colin de Verdiere [5] proved Andreev-Thurston 
Theorem in the special case of tangent circle packing, i.e., 0y = for each ij £ 
E. Chow-Luo [9] proved Andreev-Thurston Theorem of intersection angle circle 
packing, i.e., 0^ £ [0, ?] for each ij £ -E. Variational principles for polyhedral 
surfaces including the topic of circle packing are studied systematically in Luo [16j . 
Many energy functions are derived from the cosine law and its derivative. Tangent 
circle packing is generalized to tangent circle packing with a family of discrete 
curvature. For exposition of this work, see also Dai-Gu-Luo [10] . We follow this 
variational approach to prove Theorem [2] of inversive distance circle packing, i.e., 
Iij £ [0, oo ) for each ij £ E. 

2. Circle packing in Euclidean geometry 

2.1. Space of radius vectors. As stated in Introduction, to obtain a Euclidean 
cone metric, we need conditions on a radius vector. In this section we investigate 
the space of radius vectors for one triangle in Euclidean geometry. Consider a 
topological triangle Aijk £ F. Assign a number l^j to the edge ij £ E, where Uj is 
the distance of the centers of two circles with radii r$ , Tj and inversive distance ly . 
Therefore is a function of rj, rj. In fact, in Euclidean geometry 

Uj = ytrf +r| + 2r t r 3 I l3 . 

Similarly, we obtain Ijk , hi ■ 

Lemma 3. The set lZ E k := { (r^ , rj , ) £ R^q | lij,ljk>lki satisfy triangle inequal- 
ity} is a connected and simply connected open subset in R>q. 

Proof. Let pe '■ R>o — > R> be the map sending (n, rj, r k ) to (lij, Ijk, Iki)- We claim 
that Pe is a smooth embedding. First the determinant of the Jacobian matrix of 
Pe is nonzero. Next, we will show pe is one-to-one. 

That means we need to show if the following equations has a solution (rj, rj, r k ), 
then the solution is unique. 

r* + rl + 2I jk r jrk = l%, 

rl + rf + 2I hl r k r l = l 2 kl , 
r 2 l +r) + 2L l] r l r,=l%. 

Solving the first equation for rj, we get rj = —Iju^k ± \Jij]h ~ l) r fc + Since 
Ijk > and rj > 0, we must have rj = —Ijkrk + \J {1% ~ l) r it + Hence rj is 

uniquely determined by r k . By the same argument r, = — ifej^fe + \J (Jm ~ l) r ?J + I'tv 
After substituting n,rj in the third equation, we need to show it has a unique 
solution for r&. Consider the function 

f(t) = rf(t) + r 2 j(t) + 2I m (t)r 3 (t) - 1% 

which is obtained by replacing rk by t for simplicity of notations. Now 

f{t) = (2r t + 21 l3 rj)r' l (t) + (2r 3 + 2i ii r i )^(t). 
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And we have 

m = -h ^ 



If Iki S [0, 1], then each term in r[(t) is non-positive. Thus r[(t) < 0. 
If hi S (1, oo), then 

By the same argument, we have r^-(t) < 0. Since 2^+2/^^ > and 2^+2/^^ > 
0, then /'(i) < holds. Therefore f(t) = has a unique solution. 

To determine the image pb(R5, ), we consider the image of the boundary of R> . 
If j-j = 0, then r/j = hi, fj = lij. Thus 

ljk = r] +4 + 2I jkrj r k = 1% + & + -ll,,!,,!,,. 

On the other hand, if rt > 0, then r k < hi, rj < kj. We have 

Therefore the image /9b(R> ) is the subset of vectors (lij, ljk, hi) satisfying 

tfk < lij + Iki + 2Ijkkjhi, 

l>li < + 'ij + Zhiljkhj, 

'ij < Iki + 'jfc + ^hjhiljk- 

In fact, ^^(M^q) is a cone bounded by three surfaces. In R> , the intersection of 
any of the two surfaces is a straight line. The three lines of intersection are 

kj = and ljk = hi, 

Ijk = and Iki = kj, 

hi — and kj = ljk- 
Let C be the subset of R> formed by vectors (kj , Ijk, hi) satisfying the triangle 
inequality Therefore £ is a cone bounded by three planes. The three lines of 
intersection of the three planes are the same as that of the boundary surfaces of 

We claim that p^(R>o) n £ is a cone. Therefore it is connected and simply 
connected. 

In fact, if 1^ e [0, 1], then, in R> , we have 

{{lij, ljk, hi) | I>ij ^ hi ^ jk ^^ijhiljk} ^= {{hj , ljk, hi) | hj ^ hi ljk\- 

If lij 6 (1, oo), then, in M> , we have 

{(lij, Ijk, hi) | lij ^ hi hk ~t~ ^^ijhiljk} ^ {(Jij i Ijk, hi) | lij ^ hi ~t~ ljk\- 

To determine the shape of /^(K^q) n C, there are four cases to consider. 
If lij, lj k , hi G (l-oo), then pe^-Io) H £ = C 

If lij E [0, 1], ljk, hi 6 (l-oo), then p^R?^) n£ is a cone bounded by one surface 
and two planes. 

If hj, ljk 6 [0, 1], hi S (l.oo), then pb(R>o) ( ~ i >C is a cone bounded by two surfaces 
and one plane. 

ttlij,ljk,hi G [0,1], then p £ (M^ )n/:-pB(R3 ). 
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To finish the proof of the lemma, by definition, lZ E k = p b 1 (pb(R5 >0 ) n£). Since 

□ 



Pe is a smooth embedding, lZ E k is connected and simply connected 



In the construction of the energy function, we will use a natural variable u = 
(ui, ^2, u\v\) where Ui = m(rj). Let te '■ R 3 — > R>o be the map sending 
(iti, Uj, Uk) to (ri,rj,rk). Since te is a diffcomorphism, we have 

Corollary 4. T/ie setU 1 ^ := T^ 1 (72.^ fc ) is a connected and simply connected open 
subset in R 3 . 

To prove Theorem [2] we construct a concave energy function on U % E for each 
triangle Aijk £ F. We obtain the local rigidity of inversive distance circle packing 
from this energy function. We can not prove the global rigidity due to the fact 
that U l ^ k is not always convex in this general setting. In fact, kj, Ijk, hi satisfy the 
triangle inequality is equivalent to 

(1) (Zjj + Ijk + hi)(ljk + hi — hj)(ljk + hj — hi)(hj + Ijk — hi) > 0. 

Substituting the functions hj,ljk,hi in terms of Ui,Uj,Uk into inequality ([T]), we 

he in 

1 -1k , i-4 , 



see that U^ k is defined by the inequality 



(2) 



2(ljklki + hj) , 2(ifcjXy + Ijk) | 2(lijljk + Iki) q 



g«i+Mj gtlj+Ufc gMfc+Mi 

For example, when 7y = 7j& = 2^ = 2, the inequality ([2]) determines a subset 
which is not convex. 

In Andreev- Thurston Theorem, the space of cone angle is determined by con- 
tinuity argument. Since we can not prove the global rigidity, we are not able to 
determine the space of cone angle in the case of inversive distance circle packing. 

2.2. Energy function. To prove Theorem[2j we will construct an energy function 
on the space of radius vectors. First, we construct an energy function on the space 
for one triangle. To simplify notations, consider the triangle A123 £ F and the 
map 

- £ R 3 >0 

(wi, 1*2,1/3) i-> (ri,r 2) r 3 ) h-> (I1J2J3) >-> (ai,a 2 ,a 3 ) 
where rj = e Ui for i = 1, 2, 3. For three fixed positive numbers J12, 123, 131, we have 
If = r 2 j + r\ + 2IjkrjTk for {i, j, k} = {1, 2, 3}. For simplicity of notations, we set 

(h, h, h) — ('23, hi, ^12)- And a\, «2, 0:3 are inner angles of a Euclidean triangle with 

— 1? -\-i 2 

edge lengths h,h,h- By the cosine law, cosa^ = — ~ f° r {hj^} — {1, 2, 3}. 
Therefore otx, «2, cn 3 are functions in terms of iti, 1*2, U3. 

Lemma 5. The Jacobian matrix of functions ai,a2,Q;3 in terms o/iti,U2,U3 is 
symmetric. 

Proof. Due to the cosine law of a Euclidean triangles, «i,a2,a3 are functions of 
£1, 12, h- The following formula can be proved by direct calculation. For {i, j, k} — 
{1,2,3}, 
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dai 
da 2 
da 3 




dh 
dh 
dh 



Next, by the construction l\ — r| + r\ + 21j k TjT k , for {i, j, k} = {1, 2, 3}, where 
Ijk is a constant. Differentiating the two sides of the equality, we obtained 



dh = 



U + r 2 



2l t r 3 



-drj + 



2kr k 



-dr k . 



Finally, r,; = e Ui implies Tidui = dri, for i = 1, 2, 3. 
Combining the three relations, we have, for {i, j, k} 

da i 
da 2 
da 3 



{1,2,3}, 




n 







1 du\ 









du 2 








rz 1 


V du 3 



2l 3 n 



We write the above formula as 




To show the Jacobian matrix of functions ai,a 2 ,a 3 in terms of ui,u 2 ,u 3 is 
symmetric is equivalent to show the matrix TV is symmetric. 
By the cosine law, we have 



AN 



-211 

q + q- 

q + q- 



I 2 



q + 1? 



q + q 




q + ri-4 



To simplify notations, let a :~ lf,b :— q,c :— q,x :— 
Note that ax + by + cz = 0. Thus 



1 

¥ 

H 








i 





- 






i- 


rl 











i- 


r\ 




2 2 




r 2~ r 3 



o \ 




/ 



q+ri 



, y ■= ■ 



-,z := 




—2a a + b — c c + a — b 
! V= I a + b-c -2b b + c-a 
c + a — b b + c — a —2c 

2a— (a+6— c)y+(c+a— b)z c—a—b—2ax — (c+a—b)z b—a — c+2ax+(a+b—c)y 
c—a — b+2by+(b+c—a)z 2b—(b+c—a)z+(a+b—c)x a — b—c—2by — (a+b—c)x 
b—c—a—2cz—(b+c—a)y a — b — c+2cz+(c+a — b)x 2c—(c-\-a—b)x+(b+c—a)y 

This matrix is symmetric due to the equality ax + by + cz = 0. 



□ 
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Lemma 6. The Jacobian matrix of functions a\,OL2,Ois in terms of ui, 112,113 has 
one zero eigenvalue with associated eigenvector (1, 1, 1) and two negative eigenval- 
ues. 

Proof. Using the notations in Lemma [5j we need to show that the matrix AN has 
one zero eigenvalue and two positive eigenvalues. Since 47V is symmetric and the 
sum of entries in each row is zero, we can write AN as 





where B\, B 2l B3 are functions of a, b, c, x, y, z. 
The characteristic equation of AN is 

A(A 2 + 2{B 1 +B 2 + B 3 )X + 3{B ± B 2 + B 2 B 3 + B3B1)) = 0. 

There is a zero eigenvalue and the associated eigenvector is (1, 1, 1). We claim that 
B\ + B 2 + B3 < and BiB 2 + B 2 B 3 + B3B1 > 0. These two inequalities imply 
that the matrix AN has two positive eigenvalues. In the following we verify the two 
inequalities. 
First, 

B v + B 2 + B3 < 

trace of AN = -2{B X + B 2 + B 3 ) > 

a + b + c + (b - c)x + (c - a)y + (a - b)z > 

<=> l{ + il + it + {il - ilf-^ + {il - + Ql - ©V > 

l l '2 L 3 

Zj + l 2 + Z 3 

l\l\ + — l\ — l\ 2 ^1 + ffil ~ ^1 ~ ^3 2 1 ffil + ffil ~ ^1 ~ ^2 , 

+ ;2;2 r l 72/2 T 2 + ;2;2 ' 

t 2 t 3 L 1 L 3 'l l 2 
; 2 , /2 , i2 , c: „2 , n „2 , rj „2 



2 



*i + ^2 + % + Sirl + S 2 r^ + S 3 r| > 0. 



For {i, j, fe} = {1, 2, 3}, since i? = r| + r 2 . + 2Ij k rjr k and 7^ e [0, oo), we have 
Z; > rj . 

Without loss of generality, we may assume l\ > l 2 > I3. Therefore Si > and 
S 3 < 0. 

Case 1. If S 2 > 0, then 

l\ + l 2 2 + ll + Sir? + S 2 r 2 + S 3 r 2 
>l 2 + l 2 + l 2 + S 3 r 2 

2 1 ;2 1 ;2 , 77 ;2 



> + h + ^3 + "3^2 



72/2 , 72/2 _ 74 _ /4 
_ ;2 1 ;2 1 ;2 , t .3 t l ' l 3 l 2 '1 l 2 ;2 
— l l T '2 T '3 "I 7272 f ' 

n'2 

1 



2/2 '2 
l'2 

2 ]2\i2 1 07272 1 72/2\ 



2 (Gf-^2+2^+^3)>0 



due to the assumption Zi > Za- 
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Case 2. If S 2 < 0, then 

l\ + l\ + l\ + E ir f + S 2 r 2 + S 3 r 2 

> l\ + l\ + l\ + E 2 r 2 2 + S 3 r 2 

> l\ + l\ + l\ + S 2 / 2 + E 3 l\ 




+ — — 2 + W% ~ ^1 ~ ^2 

,0 ,0 ^3 "T 




= 72^1 + h + J 3 )(/i + h - h)(h + h- h)(h + h-h)>Q 



i 



due to the triangle inequality. 
Thus B 1 +B 2 + B 3 < holds. 

To prove B\B 2 + B 2 B 3 + B 3 B\ > 0, wc substitute back 



Bi = AN 23 = a-b-c-2by-(a + b- c)x, 
B 2 = 4iV 3 i = b - c - a - 2cz - (b + c- a)y, 
B 3 = 4iVi2 = c — a — b — 2ax — (c + a — b)z 



to obtain 

B\B 2 + B 2 B 3 + B 3 B\ 
= 2ab + 2bc + 2ac - a 2 - b 2 - c 2 
+ (b 2 - c 2 - 3a 2 + 6ab + 4ac)xy 
+ (c 2 - a 2 - 3b 2 + 66c + Aab)yz 
+ (a 2 - b 2 - 3c 2 + 6ac + 46c)zx 

+ 2a(a + b- c)x 2 + 26(6 + c - a)y 2 + 2c(c + a - 6)z 2 

+ 2(a + 6 + c) (ax + by + cz) 
= (2ab + 26c + 2ac - a 2 - 6 2 - c 2 )(xy + yz + zx + 1) 

+ 2(ax + by + cz)((a + b — c)x + (6 + c — a)y + (c + a — b)z) 

+ 2(a + b + c)(ax + by + cz) 

= (2ab + 26c + 2ac - a 2 - 6 2 - c 2 )(xy + yz + zx + 1) 

due to the fact ax + by + cz = 0. 
Since 



= (h + h + h)(h + h- h)(h +h- h)(h + h- h) > 0, 

we have 

B X B 2 + B 2 B 3 + B 3 B X > 
xy + yz + zx + 1 > 



(r\ r 2 2 )(r 2 3 r\)l\ + (r 2 2 r 2 3 )(r\ r\)l 2 2 + (r 2 - r\)(r\ r\)l 2 3 + l\l 2 2 l 



2ab + 26c + 2ac - a 2 
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To verify inequality (J3|), first note, since If = rf + rf + 2IjkrjTk and Ijk & [0, oo), 
we have If >rf+rf for {i, j, k} = {1, 2, 3}. 

Since the coefficient of If in the left hand side of (|3|) satisfies 

(r? - rf)(rf rf) + iflf > (rf - rf)(rf - rf) + (rf + rf)(rf + rf) 
= 2rfrf + 2rfr\ > 0, 

by replacing If by r\ + r\ , we see that the left hand side of ([3|) is not less than 

(4) (rf - rl)(rl - r\)(rf + rf) + (rf r 3 2 )(r 2 - rf)lf + (rf - r\)(rf r%)l% 

+ irl+rl)lfll 

Since the coefficient of If in ((4]) satisfies 

(rf ri)(rf rf) + (rf + rf)lf > (rf - rf)(rf - rf) + (rf + rf)(rf + rf) 

= 2r\rf + 2rfr\ > 0, 

by replacing If by rf + rf, we see that (j4j is not less than 

(5) (rf - rf)(rf r\)(rf + rf) + (rf rf)(r\ rf)(rf + rf) 

+ (rf rf)(rf rf)lf + (rf + rf)(rf + rf)lf. 

Since the coefficient of If in ((8]) is positive, by replacing If by rf+rf, we see that 
([5]) is not less than 

(6) (rf - rf)(r\ rf)(rf + rf) + (rf rf)(rf rf)(rf + rf) 

+ (rl ~ r\)(rf - r\)(r\ + rf) + (rf + rf)(rf + rf)(rf + rf). 
Since © is equal to 8rfrfrf, we see that B X B 2 + B 2 B 3 + B-^Bx > 0. 

□ 

For the special case of intersection angle circle packing, Lemma [6] was proved in 
Chow-Luo [9](Lemma 3.1). In that case, due to Thurston [21], the monotonicity 
of angles holds i.e., < 0, ^ > for j =/= i. Combining with linear algebra, 
this property implies Lemma [6] But in the general case of inversive distance circle 
packing, the geometric picture does not hold. For example, when li = 2, l 2 = 2, l 3 — 
3, r i = r 2 = r 3 = 1, the matrix N is 

9 I _ 9 

7 I J 

\ _9 9 3 
4 4 2 

with eigenvalue 0, |. 

Since, by Corollary |4j the space U]jf 3 of vectors (ui,u 2 ,u 3 ) is connected and 
simply connected, Lemma [5] and Lemma [6] imply 

Corollary 7. The differential 1-form *^ i=1 ondui is closed. For any c 6 ^i; 23 , 
the integration w(ui, u 2 , u 3 ) = J^ Ul ' U2 ' U3 ^ ^ 3 =1 ondui is a concave function on U^f 3 
and satisfying, for i = 1,2,3, 

dw 
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Proof of Theorem^ in Euclidean geometry. Let's prove the local rigidity of inver- 
sive distance circle packing in Euclidean geometry. For a closed triangulated 
weighted surface (£, T, /) with ly S [0, oo) for each ij E E, let Ue be the open sub- 
set of KJ y l formed by the vectors u = (u\, 112, .., u\v\) satisfying (ui,Uj, Uk) E U l J, 
whenever Aijk E F. 

For a radius vector r = (n, rg, ••, T| v"|) 5 we may rescale it to get cr while the 
cone angles remain the same. Under the rescaling, u changes to u + ln(c)(l, 1, 1). 
Therefore we only consider the intersection Ue H P where P is the hyperplane 

defined by X)'=i M * = 0- 

Let a,i be the cone angle at vertex i under the Euclidean cone metric. Then 
Si=i a i — tbe sum of all inner angles = 7r|.F|. Therefore the space of possible cone 
angles is contained in the intersection R^j n Q where Q is the hyperplane defined 

by a i = n \ F \- We have a ma P £ : u e n p -> R^J n Q. 

By Corollary [71 for each triangle Aijk S -F, there is a function w(ui,Uj,Uk). 
Define a function VF : — > R by 

W(u 1 ,U 2 , ...,«|V|) = J! w ( u i, u ji u k) 
AijkeF 

where the sum is over all triangles in i* 1 . 

By ([7|. equals the sum of inner angles having vertex i, i.e., the cone angle 
at vertex i. Therefore, when restricted on Ue H P, the gradient X7W is the map 

£:U E nP^ rLTJ n q. 

We claim that the restriction of W on Ue HP is strictly concave. By the definition 
of W, we get the relation between the Hessian matrixes 

H(W)= ^2 H(w(u ilUj ,u k )), 

AijkeF 

where we think w(ui, Uj, Uf.) as a function on Ue- 
For any x E Ue H P, 

xH(W)x T = xH(w{ui,Uj,Uk))x T 

AijkeF 

= 2J {x ll x J ,x k )H(w{u l ,u j ,u k ))(x l ,x : j,Xk) T < 

since H(w(ui, Uj, Uk)) is negative semi-definite for each Aijk E F. 

Thus iff(W)i T = implies (xi,Xj,Xk)H(w(ui,Uj,Uk))(xi,Xj,Xk) T = for each 
Aijk. By Lemma |6l (xi, Xj , Xk) is an eigenvector of zero eigenvalue. Therefore 
{xi,Xj,Xk) = Cijk{^, 1, 1) for some constant c^fc. If two triangles Azjfc and Ai'j'k' 
share a vertex, then c^fc = Ci/yy. This implies that cy/s must be a constant c 
independent of triangles Aijk. Therefore x = c(l, 1, 1). Since x E P, i.e, J2 x i = 
0, we have c — 0. Thus cc = (0,0,..., 0). This shows that the Hessian H(W) is 
negative definite when restricted ouUe HP. 

The local rigidity of inversive distance circle packing in Euclidean geometry 
follows from the following lemma in analysis 

Lemma 8. If X is an open set in R n and the Hessian matrix of f is positive 
definite for all x E X , then the gradient V/ : X — > R n is a local diffeomorphism. 

□ 
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3. Circle packing in hyperbolic geometry 

3.1. Space of radius vectors. As stated in Introduction, to obtain a hyperbolic 
cone metric, we need conditions on a radius vector. In this section we investigate 
the space of radius vectors for one triangle in hyperbolic geometry. Consider a 
topological triangle Aijk G F. Assign a number Zy to the edge ij G E, where Zy is 
the distance of the centers of two circles with radii r.; , rj and inversive distance 7y . 
Therefore Zy is a function of r^, rj. In fact, in hyperbolic geometry, Zy is determined 
by coshZy = coshr^ coshr-,- + Jy sinhr^ sinhr.,-. Similarly, we obtain Ijk, Iki- 

Lemma 9. The set IZ 1 ^ :— { (r^ , rj , r& ) € R5,q | lij,ljk,lki satisfy triangle inequal- 
ity} is a connected and simply connected open subset in R> - 

Proof. Let pn ■ K>o —> R>o be the map sending (r^, rj, rfc) to (Zy, Ijk, Iki)- We claim 
that ph is a smooth embedding. First the determinant of the Jacobian matrix of 
Ph is nonzero. Next, we will show pu is one-to-one. 

That means, for fixed numbers Jy , Ijk, Iki G [0, oo) and a fixed hyperbolic trian- 
gle, there is at most one configuration of three circles centered at the vertexes of 
the triangle and having the prescribed inversive distance iy , Ijk , hi ■ 

The statement is true in Euclidean geometry as proved in the proof of Lemma [3] 
We use the result in Euclidean geometry to prove the similar result in hyperbolic 
geometry. Let's assume there are two configurations of three circles in hyperbolic 
plane with the fixed triangle and inversive distance. As in Figure [2j the stereo- 
graphic projection sends the two configurations to Euclidean plane C. Since the 
stereographic projection preserve the inversive distance, there will be two differ- 
ent configuration of three circles with the same Euclidean triangle and the same 
inversive distance. This is impossible. 

To determine the image /^(R^q)) we consider the image of the boundary of R> . 
If n = 0, then r fe = lki,Tj = Zy . Thus 

cosh Ijk — cosh rj cosh r^ + Ijk sinh rj sinh r^ = cosh Zy cosh Iki + Ijk sinh Zy sinh 1^ ■ 

On the other hand, if > 0, then < Iki, fj < Zy . We have 

coshljk < coshZy cosh/fci + Ijk sinhZy sinhZy. 

Therefore the image pn (R>o) is the subset of vectors {Uj, Ijk, Iki) satisfying 

coshZjfe < coshZy coshZfei + Ijk sinhZy sinhZy, 

cosh 1^ < cosh Ijk cosh Zy + Iki sinh Ijk sinh Zy , 

cosh 1^ < cosh Iki cosh Ijk + hj sinh Iki sinh Ijk ■ 

In fact, ^^(M^q) is a cone bounded by three surfaces. In R> , the intersection of 
any of the two surfaces is a straight line. For example, if 

cosh Ijk = cosh 1^ cosh Iki + Ijk sinh Zy sinh Iki 

and 

cosh 1^ — cosh Ijk cosh Zy + Iki sinh Ijk sinh Zy , 
then the sum of the two equation gives 
= (coshlki + cosh Ijk) (cosh Zy — 1) + (7jfc sinhZy + 1^ sinhZjfc) sinhZy 

= (coshlki + coshZjfe)2 sinh 2 + (/jfeSinhZ^i + 7^ sinh Ijk) 2 cosh sinh 
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The only possibility is sinh 4f = 0. Therefore kj = and ljk = hi- 
Hence the three lines of intersection are 

lij = and Ijk = hi, 

Ijk = and l ki = 1^, 

hi = and l tj = l jk . 
Recall that £ is the subset of K.5.Q formed by vectors (kj, ljk, hi) satisfying the 
triangle inequality. Therefore £ is a cone bounded by three planes. The three lines 
of intersection of the three planes are the same as that of the boundary surfaces of 
Pi/(K 3 >o)- 

We claim that pn (R>o) H £ is a cone. 

In fact, if lij £ [0, 1], then, in R> , we have 

{{hj,ljk,hi) | coshkj < coshZ fc i coshL, fe + ly sinhZ fci sinhZ ife } 

— ^jk, hi) I lij hi ~l~ Ijk}- 

If lij £ (1, oo), then, in K> , we have 

{(lij, ljk, hi) I cosh/jj < cosily cosh^ fe + Fj sinhZ fc! ; sinhZ^} 

3 {(hj, ljk, hi) | ^ij < hi + ljk}- 

Applying the same argument in Euclidean geometry, we see ph(^> ) fl £ is a 
cone. Therefore IZ 1 ^ = p^ 1 (ph(M.> ) n £) is connected and simply connected. 

□ 

In the construction of the energy function, we will use a natural variable u = 
(ui, u 2 , u\v\) where Ui — lntanh^-- Let tjj '■ R 3 — > R>o be the map sending 
(ui,Uj,Uk) to (ri,rj,rk). Since th is a diffeomorphism, we have 

Corollary 10. The set lA l ^ k :— r^ 1 (IZ 1 ^ ') is a connected and simply connected 
open subset in R2.Q. 

To prove Theorem [21 we construct a strictly concave energy function on for 
each triangle Aijk £ F. As in the case of Euclidean geometry, the space U l ^ k is 
not always convex (for example, when = Ijk = hi = 2), we can not obtain the 
global rigidity. Therefore we are not able to determine the space of cone angles in 
the case of inversive distance circle packing. 

3.2. Energy function. To prove Theorem [51 we will construct an energy function 
on the space of radius vectors. First, we construct an energy function on the space 
for one triangle. To simplify notations, consider the triangle A123 £ F and the 
map 

Ujj 23 -> n 1 ™ - £ -» R 3 >0 
(ui,u 2 ,u z ) h-> (ri,r 2 ,r 3 ) i-> (h,h,h) >-> (ai,a 2 ,a 3 ) 

where r% = \zrw t for i = 1,2,3. For three fixed positive numbers I\ 2 , I 23 , 
we have cosh^ — cosh r.; cosh + sinh rj sinh for fc} = {1,2,3}. For 
simplicity of notations, we set (h,l 2 ,h) — (l 2 3,hi,h 2 )- And a\,a 2 ,a 3 are inner 
angles of a hyperbolic triangle with edge lengths l\,l 2 ,l 3 . By the cosine law, 

— cosh li + cosh lj cosh 
sinh L- sinh Zfc 
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for {i,j, k} = {1, 2, 3}. Therefore a\, 012, 03 are functions in terms of u\, U2, 113. 

Lemma 11. The Jacobian matrix of functions a±, 0:2,0:3 in terms of U\, 112,113 is 
symmetric. 

Proof. Due to the cosine law of a hyperbolic triangles, 01,0:2,0:3 are functions of 
h,l2, h- The following formula can be proved by direct calculation. For k} = 
{1,2,3}, 




-1 



sina-i sinh/, sinh/fc 




Next, by the construction cosh ^ = cosh rj cosh rk+Ijk sinh rj sinh r^, for {i,j, k} = 
{1,2,3}, where I{j is a constant. Differentiating the two sides of the equality, we 
obtained 

„ — cosh r-fc + cosh k cosh r,- , — cosh r, + cosh U cosh r^ , 

dl i = • u 7 r ~C dr J H • u 1 ■ Z drk - 

sinn li sinn r, sinn t, sinn r^ 

Finally, r, = In y^fr^ implies sinhridui = dri, i = 1, 2, 3. 
Combining the three relations, we have, for {i, j, £;} = {1, 2, 3}, 




-1 



sin ai sinh Zj sinh Ik 




-1 

COS «3 

cos a.2 cos «i 
cosh/i coshr2 




— cosh r3 + cosh h cosh n 

sinh I2 sinhri 

— coshr2 + coshl 3 coshri 

sinh/3 sinhri 



sinhZi sinhr2 


cosh n + cosh I3 cosh r 2 



— cosh r2 + cosh /1 cosh r3 \ 

sinhJi sinh T3 

— coshri + cosh/ 2 coshr 3 

sinh I2 sinh r3 





We write the above formula as 
da 1 
da 2 
da 3 



-1 



sin Oi sinh L sinh Ik 




where M is a product of four matrixes. 

To show the Jacobian matrix of functions 011,0.2,013 in terms of 111,112,113 is 
symmetric is equivalent to show the matrix M is symmetric. 
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To simplify notations, let a :— cosh I i,b := cosh^c := cosh/3,2: := coshri,y := 
coshr2, z := coshr3. By the cosine law, we have 

1 — a 2 ab — c ca — b \ / a i_ 1 \ / ay — z az — y 
M = I ab - c 1 - b 2 be - a \ \ -J-^ 6a; - z bz - x 



ca 



b be — a 1 — 




7^1 

We check that M is symmetric, for example, by showing 

ac — b be — a 
M12 = M 2 i = z 3 -x = -y. 

c z — 1 c z — 1 




□ 



Lemma 12. TTie Jacobian matrix of functions a±, 02,0(3 in terms of Ui,u 2 ,U3 is 
negative definite. 

Proof. Using the notations of Lemma [Til we need to verify that the matrix M is 
positive definite. The determinant of M is positive since the determinant of each 
factor is positive. Note that the second matrix is a Gram matrix up to the negative 
sign. Since, by Corollary QJJ the set Uj^ 3 is connected, to show M is positive 
definite, it is enough to check M is positive definite at one vector in Uj^ 3 . 

It is equivalent to pick up a vector (r±,r2,r3) such that the resulting li,h,h 
satisfy the triangle inequality. We pick up (ri,r2,r3) = (s,s,s). In the following, 
we show that such kind of number s exists. 

For fixed I12 , 123 , 131 € [0, 00), we have cosh ^ = cosh rj cosh rk+Ijk sinh rj sinh r^ = 

cosh 2 s+Ijk sinh 2 s > cosh 2 s, for k} = {1, 2, 3}. Therefore sinh li > \J cosh 4 s — 1 
for i = 1,2,3. 

That l\, I2, 13 satisfy the triangle inequality is equivalent to 



cosh(Zj + Ij) > coshZfe 



■ cosh h cosh/. 



sinh U sinh Ij > cosh Ik — cosh s + 



Iij sinh s 



- cosh 



cosh 4 s — 1 > cosh 2 s + I;, sinh 2 , 



Let a; := cosh 2 s. The above inequality is equivalent to (2x — 1^ + l)(x — 1) > 0. 
Since x > 1, we need x > Jl3 2 1 . 

Therefore £1, 12, 13 as functions of (s, s, s) satisfy the triangle inequality if cosh 2 s > 
max{-^-, ^s_zl, ^i}. 

Now it is enough to check that M is positive definite at such a vector (s, s, s) 
when s is sufficiently large. In this case, 



M 



sinh/i 






sinh I2 
sinh £3 



— 1 cos «3 cos «2 
cosa3 —1 cosai 
cosa2 cosai —1 



— 1 + coshZi — l + cosh^i \ 

sinh/i sinh/i 
— l + cosh^2 — l + cosh^2 











sinh l 2 

— 1 + cosh I3 — 1 + cosh I3 



sinh h 



sinh is 



sinh/2 




coshs := Mi cosh s. 
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M is positive definite is equivalent to that M\ is positive definite. Since M\ only 
involves h, I2, h, that Mi is positive definite is a property of a hyperbolic triangle. 
In fat, that M\ is positive definite was proved in the paper Guo-Luo [12) (Lemma 
4.4), by showing that its leading principal lxl and 2x2 minor are positive. For 
completeness, we include a proof here. 

That the leading principal lxl minor is positive is equivalent to 

— l + cosh/ 2 — l + cosh/ 3 

cos a 3 — — 1- cos a 2 r-r-, > 0. 

smii 1 2 Sinn 43 

Replacing a's by I's by the cosine law, we obtain 

. . , , cosh 2 I2 + cosh 2 I3 + cosh I2 + cosh I3 

(8) cosh/i > — — — — — — . 

2 cosh <2 cosh t 3 + cosh 12 + cosh {3 

Since li > \l 2 — Z 3 |, therefore cosh^ > cosh(/ 2 — h)- To show © holds, it is 
enough to check 

fn . . . . . ^ cos h 2 h + cosh 2 l 3 + cosh h + cosh l 3 

(9) cosh l 2 cosh h - smh h smh l 3 > — — — — ■ — — ■ — — . 

2 cosh 62 cosh 43 + cosh 12 + cosh 1% 

We simplify the notations by introducing a := cosh/2 > 1, b := cosh/3 > 1. Then 
© is rewritten as 



(10) ab-^ia 2 - 1)(6 2 - 1) > 



2ab + a + b 
TOf is equivalent to 

a 2 + b 2 + a + b 



(11) ab- % ,V "T > y/(^-lW-l). 

2ab + a + b 

Since a > 1, b > 1, the left hand side of (fTT|) is positive. To show (fTl) holds, we 
square the two sides and simplify. We have 

(ab 4 + a% - a 3 b 2 - a 2 b 3 ) + (a 4 + b 4 - 2a 2 b 2 ) + (a 3 + b 3 - ab 2 - a 2 b) > 
^=>(ab + l)(a 3 + b 3 - a 2 b - ab 2 ) + (a 2 - b 2 ) 2 > 
^(ab + l)(a + b)(a - b) 2 + {a 2 - b 2 ) 2 > 0. 

This shows that the leading principal lxl minor is positive. 

To show the leading principal 2x2 minor of Mi is positive, for simplifying 
notations, let U := "'t™ 1 *^ . Up to a diagonal matrix, Mi is 

— 1 COS Q!3 COS CK2 

cosa3 —1 cosai 

COSQ!2 COS Oil —1 
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That the leading principal 2x2 minor of M\ is positive is equivalent to 

(t 2 COSQ!3 + £3 COS«2)(^l COS 013 + t 3 cosai) — (— t 2 + £3 cosai)(— 1\ + t 3 cos a 2 ) > 
(cos 2 a 3 — l)tit 2 + (cos ot\ cos a 3 + cos a2)t2t 3 + (cos a 2 cos a 3 + cos a\)tit 3 > 
■<=>■ — sin 2 a 3 t\t 2 + cosh l 2 sin a\ sin a 3 £ 2 £3 + cosh l\ sin a 2 sin a 3 t\t 3 > 
■<=> cosh ^2 sinh l\t 2 t 3 + cosh Zi sinh l 2 tit 3 > sinh l 3 t\t 2 
^ ^ cosh ^2 sinh ?i ^ cosh Zi sinh l 2 sinh/3 

ti t 2 t 3 

<==^ cosh/2(l + cosh^i) + coshZi(l + COSI1Z2) > 1 + COSI1Z3 
■<=>■ (cosh li + cosh I2 — 1) + (2 cosh Zi cosh Z2 — cosh l 3 ) > 0. 

The cosine law and the sine law of a hyperbolic triangle are used in above calcula- 
tion. 

The last inequality is true since coshJi + cosh/2 > 1 and 2 cosh li cosh l 2 > 
cosh li cosh l 2 + sinh li sinh l 2 = cosh(/i + l 2 ) > cosh l 3 . 

□ 

For the special case of intersection angle circle packing, Lemma Q2] was proved 
in Chow-Luo [9] (Lemma 3.1) by applying the geometric picture of monotonicity of 
angles (an observation due to Thurston [21]). But in the general case of inversive 
distance circle packing, the geometric picture does not hold. For example, when 
li =1,1-}, = 2, l 3 = 3, ri = r 2 = r 3 = 1, the matrix M is 

6.08 0.49 -2.94 
0.49 6.08 -2.94 
-2.94 -2.94 22.11 

with eigenvalue 23.15,5.59,5.53. 

Since, by Corollary [T0l the space U^ 3 of vectors (ui,u 2 ,u 3 ) is connected and 
simply connected, Lemma [Til and Lemma [121 imply 

Corollary 13. The differential l-form X^=i otidui is closed. For any c E Uj^ 3 ', 
the integration w(ui, u 2 , u 3 ) — J^" 1,tl2 '" 3 ) J^? =1 otidui is a strictly concave function 
on Uj^ 3 and satisfying, for i = 1,2, 3, 

(12) d^r ai - 

Proof of Theorem^ in hyperbolic geometry. Let's prove the local rigidity of inver- 
sive distance circle packing in hyperbolic geometry. For a closed triangulated 
weighted surface (S, T, I) with € [0, 00) for each ij E E, let Uh be the open sub- 
set of Rl^l formed by the vectors u — (u\, u 2 , .., U\y\) satisfying (ui,Uj,Uk) E Ujj 
whenever Aijk € F. 

By Corollary 1131 for each triangle Aijk E F, there is a function w(ui,Uj,Uk). 
Define a function W : Uh — * K by 

W(u 1 ,U 2 , ... ,U\ V \) = U>(Ui,V,j,Uk) 
AijfcGF 

where the sum is over all triangles in F. By Corollary [Jj2 W is strictly concave on 
Uh so that equals the sum of inner angles having vertex i, i.e., the cone angle 
at vertex i. That means the gradient of W is exactly the map sending a vector u 
to its cone angles. 
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The local rigidity of invcrsive distance circle packing in hyperbolic geometry 
holds due to Lemma [8] 

□ 
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